Importance of exchange-anisotropy and super exchange for the 
spin-state transitions in LnCo0 3 (Ln=La,Y,RE) cobaltates 
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Spin-state transitions are the hallmark of rare-earth cobaltates. In order to understand them, it 
is essential to identify all relevant parameters which shift the energy balance between spin states, 
and determine their trends. We find that A, the e g -t2 g crystal-field splitting, increases by ~ 250 
meV when increasing pressure to 8 GPa and by about 150 meV when cooling from 1000 K to 5 K. 
It changes, however, by less than 100 meV when La is substituted with another rare earth. Also 
the Hund's rule coupling J& v& is about the same in systems with very different spin-state transition 
temperature, like LaCoC>3 and EUC0O3. Consequently, in addition to A and J7a,vg, the Coulomb- 
exchange anisotropy Aj7av g and the super-exchange energy-gain AEse play a crucial role, and are 
comparable with spin-state dependent relaxation effects due to covalency. We show that in the 
LnCoOs series, with Ln=Y or a rare earth (RE), super-exchange progressively stabilizes a low-spin 
ground state as the Ln 3+ ionic radius decreases. We give a simple model to describe spin-state 
transitions and show that, at low temperature, the formation of isolated high-spin/low-spin pairs is 
favored, while in the high-temperature phase, the most likely homogeneous state is high-spin, rather 
than intermediate spin. An orbital-selective Mott state could be a fingerprint of such a state. 

PACS numbers: 71.27.+a, 71.30.+h, 72.80.Ga,71.10.-w, 71.20.-b, 75.40.Mg 



I. INTRODUCTION 



The nature of the spin-state transitions in LaCoOa 
(Fig. [1]) is the subject of controversy for decades^— 
LaCoO-3 undergoes two electronic crossovers, at Tss ~ 
50 - 100 K and Tim ~ 500 - 600 K; the first is commonly 
ascribed to a change in spin-state, while the second to an 
insulator to metal transition* 2 ^ There is a general agree- 
ment that the ground state of LaCoOa is insulating and 
non-magnetic, with Co in the low-spin (LS) ifg state. 
The core of the debate is whether the t^e^ high-spin 
(HS) stated- (5 = 2), or the t\ g e g intermediate-spin (IS) 
stated— (5=1) is thermally excited right above Tss, an d 
if the spin-state crossover involves two or more steps be- 
tween Tss and Tim- Various mixed phase a 16 ' 17 and spin- 
state superlatticesi£~— have been suggested as final or in- 
termediate step. For the high-temperature phase, a pure 
IS state has been proposed based on LDA+?7 results^ 7 - 
or phcnomenological thermodynamic analysis* 2 ^ To date, 
the controversy remains open. 

Spin-state transitions have also been reported in 
LaCoC>3 under pressure; 2 ^— as well as in several 
L11C0O3 perovskites, where Ln is Y or a rare earth (RE) 
heavier than LaSZr— In the LnCo03 series, the spin- 
state and, to a smaller extent, the insulator to metal 
transition, shifts to higher temperatures with decreasing 
RE 3+ ionic radius^Zr— Since HS and IS states are Jahn- 
Teller (JT) active, one might expect co-operative JT 
distortions in correspondence with the spin-state tran- 
sition. Remarkably, early neutrons diffraction data* 3 - 1 . - — 
indicate that the structure of LaCoOs is rhombohedral 
(space group i?3c) with no co-operative JT distortion at 
all ^i? 1- ?? while more recently a co-operative JT distor- 




FIG. 1: (Color on-line) Rhomobohedral structure of LaCo03. 
We define the x, y, z pseudo-cubic axes as those connecting 
the Co atoms; z is the vertical axis, y the horizontal axis, and 
x is perpendicular to the plane of the figure. 



tion (space group 12/ a) has been found via single-crystal 
X-ray diffraction] 14 ' 34 The presence of local JT distor- 
tions is as controversial^^ On the other hand, LnCoO*3 
perovskites with Ln^La are orthorhombic (space group 
Pbnm) and do exhibit a small, weakly temperature de- 
pendent, JT distortions 7 .^^ 

The LS— >IS scenario^ was popularized by total-energy 
LDA+J7 calculations,™ which show that for rhombohc- 
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dral LaCo03 the homogeneous IS state is almost degen- 
erate with the LS ground state, and sizably lower in en- 
ergy than the HS state . 12 ' 17 " — The LS— >IS scenario has 
been used to interpret experiments pointing to a spin 
triplet j 11 ' 14 ' 28 However, a triplet can also arise from the 
splitting of the HS state via spin-orbit interaction . 13 i 36 i 37 
Furthermore, it has been shown with LDA+J7 and unre- 
stricted Hartree-Fock that some inhomogeneous LS/HS 
phases and/or super lattices ^r— are more stable than 
the homogeneous IS state. Experimentally, neither the 
HS+LS nor other superlattices have been reported so far; 
instead, there are strong indications that X-ray absorp- 
tion spectroscopy (XAS) and magnetic circular dichro- 
ism (MCD) data 3 - as well as inelastic neutron scattering 
experiments 4 - are compatible with a mixed HS/LS phase. 
It has been pointed out that the lattice could play a cru- 
cial role in such mixed phases, by expanding around HS 
ions due to covalency effects 1 "" 3 -" 2 ^ 

While the experimental evidence in favor of LS — > HS 
crossover at T$s, perhaps with disorder phenomena, is 
growing; 3 ' 4 ' 11 the actual parameters which tip the energy 
balance in favor of a given spin-state in the real mate- 
rials have not been fully identified, and, to the best of 
our knowledge no systematic attempt to determine their 
evolution in the LnCoOs series exist. Furthermore, the 
interplay between spin-state and insulator to metal tran- 
sition, even in the simplest homogeneous scenarios, and 
the nature of the high-temperature metallic phase, are 
not fully understood. 

In the present work, we study the electronic structure 
trends of the L11C0O3 family and identify the material- 
dependent parameters responsible for the spin-state tran- 
sitions. We show that, for realistic Coulomb and ex- 
change parameters, an IS homogeneous phase is very un- 
likely in any of the considered materials, both at low 
and high temperature. We show that, besides the in- 
terplay of the e g -t2 g crystal-field A and the Coulomb 
exchange jTavg, the exchange anisotropy A^avg and the 
super-exchange energy gain Ase stabilize the LS ground 
state. Our results support a scenario in which HS-LS 
pairs, the formation of which is favored by both super- 
exchange interaction and covalency effects, arc excited 
at Tss- We find that, for all known experimental struc- 
tures, the Jahn- Teller crystal-field splitting is weak. In 
the high-tcmpcraturc regime, we find, within a homoge- 
neous scenario, a HS state; in such a state, the insulator- 
to-metal transition depends mostly on the ti g degrees of 
freedom. Our results could be used to distinguish a pure 
high-temperature HS state from a mixed HS/LS state. 



II. METHOD 

We calculate the electronic structure using two ab- 
initio approaches based on density-functional theory in 
the local-density approximation (LDA). The first is the 
Linear Augmented Plane Wave (LAPW) method as im- 
plemented in the Wien2k codej^ 8 . to obtain hopping in- 



tegrals and crystal-field splitting we use maximally lo- 
calized Wannier functions.— The second approach is the 
downfolding technique based on the Ath- Order Muffin- 
Tin Orbital (NMTO) method in the form discussed in 
Refs. By means of these two methods, we construct 
the generalized Hubbard model for the Co d bands 

H = - V t M ' ,c f c, 

/ _j rn,m' truer i'm'cr 
im^i'm'a 

+ n U™pm'p' C lma C ipa fC ip'a ,C irn'(T m CO 



2(7, cr mm pp 

Here c imo . creates an electron with spin a in the Wannier 
orbital m at site i (m = xz,yz,xy,x 2 — y 2 ,3z 2 — r 2 ). 

ipi are rotationally invariant screened Coulomb in- 
tegrals. The parameters U mpm 'p' can all be expressed 
as a function of the (screened) Slater integrals F , F 2 
and F±. We adopt the common definition for the av- 
erage direct and exchange couplings, t/ avg = Fq, and 
•^avg = nC^ 2 + -^t)- A pedagogical derivation of the 
rotationally invariant Coulomb matrix can be found, 
e.g., in Ref. lU Apart from the average exchange in- 
teraction, the exchange anisotropy plays a significant 
role. It is therefore convenient to express all parame- 
ters as linear combinations of £/ avg , J7avg = f ^avgj an d 
AJavg = J avs (l - §ff)/(l + ff); the latter measures 
the anisotropy in the exchange interactions. The ex- 
change parameters for the t2 g and e g states are then 
Jt2g = iTavg + Aj7avg and J eg = ^Tavg + 3Aj7avgj respec- 
tively; the Coulomb exchange parameters between e g and 
t 2g states are J^ z 2- r 2 xz = J^ z 2_ r 2 yz = Javg — 3Aj7avg, 



J r 



JTavg 5Aj7a V g, J^ z 2_ r ^ 

J x 2 -y 2 ,y Z = Jt 2g ] the orbital-diagonal di- 
rect exchange is Uo = Ua, vg + §,7avg- For some materials, 
we calculate the screened Coulomb integrals C/ avg and 
J7avg using the constrained local-density approximation 
(cLDA) approach^ 2 - 

We solve Hamiltonian ([1]) using different approaches. 
In section III. A we present exact diagonalization re- 
sults (atomic limit). In section III.B we show re- 
sults obtained with second-order perturbation theory 
(super-exchange energy gain). In section III.C we 
present dynamical mean-field theory calculations within 
the LDA+DMFT (local-density approximation + dy- 
namical mean-field theory) approach; for the latter we 
use a weak-coupling continuous-time quantum Monte 
Carlo^ 3 - (CT-QMC) quantum-impurity solver and work 
with the full self-energy matrix in orbital space , 40 ' 44 as 
discussed in Ref. HH. We obtain the spectral-function 
matrix by means of stochastic reconstruction^ Since 
LDA+DMFT calculations for a 5-band model with full 
Coulomb vertes^ 7 - and self-energy matrix are computa- 
tionally very expensive, a massively-parallel implemen- 
tation as presented in Ref. is essential. Structural 
data for LaCoOs arc taken from Rcfs. 32.34, for YC0O3 
from Refs. and for the other LnCo03 materials^ 8 - from 
Refs. I28ll29ll49l . High-pressure data are from Ref. [H 
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FIG. 2: (Color on-line) Left: Atomic spectral-functions at 
29 K for different values of the ratio A/ Javg- Full line, dark: 
A/J a vg ~ 1.69. Dotted line: A/J avg ~ 2.60. Full line, light: 
A/Javg ~ 2.61. Coulomb parameters: C/ a vg = 3 eV, J a v g = 
0.89 eV. Center: occupation of the e a states as a function of 
A/ Javg at low and high temperature. Right: occupation of 
the e g states as a function of AJa Vg /Ja V g at low and high 



temperature; the results are obtained for 



0.89 eV, 



A/Javg ~ 2.64 and varying F 4 /F 2 in the interval [0.1,0.8], 



III. RESULTS 



A. Atomic limit, Jahn-Teller, and constraints 



Let us examine the eigenstates of (JXJ) in the atomic 

limit, i.e., in the limit , = for i ^ i' (see Fig. [5]). 
The crystal-field states, the states which diagonalize the 
on-site matrix t 



have energy e a , with e a < £ +i; 
a = 1, 2, 3 are i2g-hke states and a = 4, 5 e s -like. In the 
atomic limit, the energy of the low-spin state (t| ff ) i s 



E L ~ 15U - 30( Javg + AJ £ 



avg l 



■ E 

a£t 2g 



The states with intermediate spin (i|g e g) can be written 
as | a, /3), where a is the ti g hole orbital and /3 the e g elec- 
tron orbital. The energy of the IS, in the limit in which 
only the density-density Coulomb terms contribute, is 



Ei 



E, 



A 



LI 



3Jav g + (23 + /)A J avg 

where A LI = £4 — 23; the factor / yields the deviation 
from the average anisotropy of the exchange interaction; 
/ = —8 for \xy,x 2 — y 2 ) and cyclic xyz permutations. 
The energy of a high-spin ifg e g is 



E H 

. LH 



E L - 8 Ja VS + 30AJ avs + 2A 



LH 



where A Ltl — (e§ + £4 — £3 — £2)/2. From energy dif- 
ferences between spin configurations we can obtain con- 
straints for the parameters. Theoretical estimates^ yield 
F-2 ~ 10.64 eV and F4 ~ 6.8 eV; with these values 



avg 



0.89 eV, and AJ~ a , 



0.07 eV - 0.08J E 



avg ■. 



other estimates^ yield slightly smaller values; by us- 
ing cLDA we obtain J avg ~ 0.7 cV; furthermore, we 
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FIG. 3: (Color on-line) Crystal-field splittings as a function 
of pressure (left), temperature (center), and of Ln ion (right). 
The corresponding volumes are also given; Vo is the vol- 
ume of LaCo03 at room temperature and ambient pressure. 
A LI = e 4 - £3 (squares) , A avg = (£4 + £5)/2-(£i+£ 2 + £3)/3 
(dark circles), and A LH = (e 4 + e 5 )/2 - (£2 + e 3 )/2 (light cir- 
cles) are shown. In the cubic limit, A LI = A LH = A avg = A. 
Multiple points for the same compound are results from dif- 
ferent structural data sets. 



find that j7a Vg varies little when La is replaced by Eu in 
EuCo03, a material which apparently exhibits no spin- 
state transition below the insulator-to-metal transition^ 
Since screened Coulomb exchange parameters can only 
be obtained within given approximations, we will discuss 
our results for a range of plausible values of Javg , ranging 
from 0.5 eV to 1 cV. 

In Fig. [3] we show the calculated crystal-field splitting. 
We calculate the splittings in LaCo03 at ambient pres- 
sure and for increasing temperature) 31 ' 32 at room tem- 
perature for increasing pressure 2 ^ 26 ' 52 and in the se- 
ries of REC0O3 compounds . 27 ' 49 ' 53 ' 54 Our results show 
that A LH - 1.6 - 1.7 eV in the full series,-^ while 
A LI ~ 1.4 — 1.6 cV. From these results we can derive the 
following conclusions. The ground state is LS if these 
conditions are met 



A LH > 4J £ 



aV g 15A^7avg ^ 2.8 J/avgi 



.LI 



> 3 Javg - (23 + /)AJ avg ~ 1.16-1.8 J„ 



where for the second condition we give the range vary- 
ing / from (average) to —8 (\xy,x 2 — y 2 ))- Let us 
consider the cubic case (A LH = A LI = A avg = A). 
In this case, the first condition is the most stringent. 
Exact diagonalization of the atomic limit Hamiltonian 
with full exchange interaction and AJ aV g/J aV g ~ 0.08 
(Fig. [2]) indeed leads to a switch between a high and 
low spin ground state at slightly smaller values, A ~ 
2.6 J aV g. From these results we conclude that, to explain 
a S = ground state in the atomic limit, J avg must 
be slightly smaller than estimated with cLDA and sim- 
ilar approaches, J avg ~ 0.6 eV; alternatively, less likely, 
Javg ~ 0.7 eV but the anisotropy should be sizably larger 
than in the free atom,— AJ avg ~ 0.12J avg . In any case, 
with the crystal field splittings in Fig. [31 we are quite 
close to the HS/LS crossover. In this parameter region, 
the ~ 100-200 meV changes in crystal-field that we find 
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FIG. 4: (Color on-line) Crystal-field splitting (in eV) of e g 
(left) and t2 g (right) states at room temperature as a function 
of decreasing ionic radius. Multiple points for the same com- 
pound are results from different structural data . 27 ' 49 ' 53 ' 54 ' 58 
LaCoOs: 12 /a (diamonds) and R3c (circles). 



(Fig. [3]) with increasing temperature or pressure could 
indeed explain alone a spin-state crossover from LS to 
HS or viceversa. For a IS ground state, the following 
conditions should be met 



2A LH_ A LI > ( 5 _( 7 _ /) ^ 
.LI 



avg 

A Li < 3Jkv g - (23 + /)AJkv g ~ 1.16-1.8 J„ 



This definitely excludes a IS ground state in the cubic 
case as, in all systems considered, the two conditions are 
never satisfied at the same time. For a LS— > IS sce- 
nario, it could be however sufficient that the IS state is 
lower than the HS (first condition). In the cubic case, 
for the parameters in Fig. [3j this could happen only for 
apparently unrealistically small Coulomb exchange j7a Vg 
(~ 0.4 eV or smaller) or unrealistically large^ 6 - anisotropy 
AJavg/iJavg (~ 0.16); in these cases the IS would, how- 
ever, be quite high in energy (~1 eV or more) above 
the LS. Thus the predictions of conventional ligand-ficld 
theory 5 -^ are in agreement with our results for the atomic 
limit. 

Let us analyze the effects of crystal distortions on the 
crystal-field splittings of e g (A54 = £5 — £4) and ti g 
(A31 = £3 — 61, A21 = £2 — £1) states, and the corre- 
sponding crystal-field orbitals. In Fig. 0] we show A54, 
A31 and A21 along the scries; we find that the split- 
ting of e g states, A54, reaches a maximum of 200 meV 
in H0C0O3, and 250 meV in tetragonal 3 ^ LaCoC>3, and 
is zero in rhombohedral LaCoOs; furthermore, we find 
that A21 is small while A54 and A31 are comparable, the 
only exception being the 12/ a LaCoOs structure; finally, 
across the spin-state transition, for available structural 
data, we find either basically no change in e g and t2 g 
crystal-field splittings (LaCoOs, i?3c) or a tiny change 
of 40 meV and no significant change in occupied orbital 
(YC0O3). Even if we assume that the crystal-field split- 
tings in Fig. @] are completely ascribed to the Jahn- Teller 
effect, they are small compared to Jahn- Teller systems 
such as manganites or cuprates^ they are compara- 
ble with splittings arising from the GdFe03-type distor- 
tion (Pbnm structures) 40 ' 59 which progressively increase 
along the LnCoC>3 series. 53 ' 54 LaCoC-3 has to be discussed 



separately. While the presence of JT distortion is still 
controversial experimentally, even the monoclinic struc- 
ture with co-operative Jahn- Teller distortion reported in 
Ref. HH leads to a small splitting of e g and t2 g states 
(see Fig. [4j, although larger than for all other LnCo03 
compounds considered here. 

The splittings of e g and t2 g states might affect the en- 
ergy balance between spin states. To quantify this effect, 
we introduce the average splitting SA = (A54 + A3i)/2. 
Then, on the basis of Fig. H A LH ~ A, A LI — A — 8 A. In 
the light of presently known structural data, remarkably, 
in the full LnCo03 series, 8A is at most 150 meV, which 
is insufficient to lead to a LS-IS scenario, which would 
require 8A ~ 1-1.7 Ja. vg or larger. Such a scenario could 
perhaps start to play a role under high pressure, if 8A 
sizably increases. 

Our results show that the anisotropy AJ& V g plays an 
important role in stabilizing the low-spin ground state 
(Fig. [2]); if we neglect it (Ji = J w% = J), an approxima- 
tion often adopted, the constraint for a low spin ground 
state becomes A LH > A J . Thus, if the e g -t2 g crystal-field 
splitting A has the values of A avg in Fig. [3l the effective 
Coulomb exchange has to be as small as J ~ 0.4 eV in 
order to obtain a LS ground state. On the other hand, if 
the anisotropy is unrealistically larger than for atomic or- 
bitals, the IS state can become the first excited state. As 
a consequence, for systems close to spin-state transitions, 
differences in total energy between spin states might be 
extremely sensitive to the approximations adopted in de- 
scribing the multiplet structure ! 60 ' 61 

The value of £/ avg does not affect the energy differences 
between spin-states in the atomic limit, but is relevant 
for the insulator to metal transition and the gap, as well 
as for super-exchange. Theoretical estimates based on 
density- functional theor y 12 ' 50 ' 62 yield ~ 8 eV in LaCoOs, 
and we find similar values with constrained LDA. Ex- 
perimental estimates based on electron spectroscopy 6 ^ 
yield sizably smaller values (~ 3.5 eV); similar values 
were adopted in many worksi 17 ' 19 ' 37 A configuration- 
interaction cluster model analysis of spectroscopic data 
yields C/ avg = 5.5 eV for the low and intermediate spin- 
state and ~ 6.8 eV for the high-spin state. Since the 
value of £/ avg affects energy scales such as the size of the 
gap and super-exchange, when necessary, in the next sec- 
tions, we present results for several J7 avg . 

In conclusion, our results show that, in the atomic 
limit, for all systems analyzed, a LS— >IS scenario is un- 
likely At close inspection, LDA+C7 total-energy calcula- 
tions are in line with this conclusion ; 12 ' 17 ~ — as they in- 
dicate clearly that the stabilization of a given local mag- 
netic state (IS or HS) strongly depends on the number 
and type of surrounding magnetic neighbors, suggesting 
that co-operative, band, or at least multi-site effects play 
a crucial role; at low concentrations, in the absence of 
magnetic neighbors, the HS appears to be the favored 
magnetic statei^ Another important aspect is the evo- 
lution of the parameters with increasing temperature. 
For LaCoOs we find that the crystal-field A avg is re- 
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duced by about 150 meV when the temperature increases 
from 5 K to 1000 K. Remarkably, magnetic susceptibil- 
ity dataii indicate that the LS— > HS activation energy 
Aae, increases with temperature; similarly the analysis 
of XAS data by means of the truncated configuration- 
interaction cluster approach^ suggests a rise of about 60 
meV from 50 K to 700 KJ^ This result cannot be ex- 
plained in the atomic limit, even including the effects of 
spin-orbit interaction, AS ■ L; we calculate the spin-orbit 
coupling A for LaCo0 3 and find A ~ 54 meV, slightly 
larger than the one assumed in model calculations^ 
but too small to affect the trends on A avg . Our results 
also exclude that the rhombohedral distortions increase 
enough with temperature to overcome the decrease in 
A avg due to the increasing volume^ 



B. Super exchange 

Away from the atomic limit, super-exchange affects 
the energy balance. However, is the super-exchange en- 
ergy gain large enough to be relevant in the spin-state 
crossover at Xss? In this section we calculate the super- 
exchange energy gain for several types of Co-Co bonds: 
LS-LS, HS-HS and LS-HS pairs; we do not consider IS-LS 
pairs, whose formation at low temperature is unlikely^ 
While differences in the energy of multiplets with fixed 
number of electrons, relevant to determine the spin-state, 
arc of the order J , virtual excitations of electrons to 
neighboring sites, relevant for super-exchange, involve 
the direct Coulomb energy IA. Thus Coulomb exchange 
anisotropy A^vg has a small effect on super-exchange; 
for simplicity we neglect it and consider density-density 
terms only. For a LS-LS pair we obtain the energy gain 



= -EE 

a=l a'=i 



U -5J + A a ., 



(2) 



where A a >, a = £a' — For a HS-HS pair, assuming 
that Co ions are in a paramagnetic state, we obtain 
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FIG. 5: (Color on-line) Super-exchange energy gain per bond, 
Asb, for a LS-LS pair (empty triangles), a HS-HS pair (filled 
triangles) and a LS-HS pair (pentagons). Left: Ase as a func- 
tion of pressure. Center: Ase as a function of temperature. 
Right: Ase as a function of Ln ion. The corresponding vol- 
umes are also given; Vo is the volume of LaCo03 at room tem- 
perature and ambient pressure. Multiple points for the same 
compound are results from different structural data sets. 



state for LaCoOa for all temperature and pressures con- 
sidered. In the series LnCoOa, the super-exchange energy 
gain becomes larger for a LS ground state for ions smaller 
than Dy 3+ . For [/ avg ~ 3.5 eV the crossing happens al- 
ready for Ln=Sm, and the energy gain per LS-LS bond 
increases sizably to 90 meV around Ln=Y. 

It is at this point crucial to evaluate also the super- 
exchange energy gain associated with the formation of a 
HS-LS bond, given by 



ae£I ls 
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+ E 



5 \A,i |2 , uM' |2 

E^ 



P 



U -J 



5 'I^+IC,/ , l$J a +lff,/ 



U -7J + A a ,, f3 U -3J + A a , t/3 



3 \J,,i |2 , ui,» |2 

E^ E^ 



U - J + A c 



(3) 



where {/?} are the np degenerate ti g levels, over which 
we average. Results for material-dependent hopping inte- 
grals are shown in Fig. \5\ We find that for £/ avg ~ 5 eV or 
larger the energy gains do not change much with increas- 
ing J from to 1 eV^ Instead, for fixed J ~ 0.89 eV, 
we find a remarkable change in Ai?g|" LS reducing f7 avg 
to ~ 3.5 eV, because the denominator sizably decrease. 
In Fig. [5]we show two significant parameter ranges. For 
t^avg ~ 5.5 eV, super-exchange favors a high-spin ground 



This energy gain is relevant in many scenarios of spin- 
state crossover. We find that A£f|~ LS /A£^|~ HS is 
~ 1.5-2 for C/ avg ~ 5.5, and it increases with decreas- 
ing ZTavg. Furthermore, all super-exchange energy gains 
calculated here vary slowly with increasing temperature. 
These results show that, even neglecting lattice relax- 
ation effectSfi it is energetically favorable to form isolated 
HS ions rather than HS-HS bonds. 

The super-exchange energy difference between a HS- 
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FIG. 6: (Color on-line) Magnetic coupling Jse for a HS-HS 
pair, as a function of pressure (left), temperature (center), and 
of Ln ion (right; empty symbols: (001) direction; filled sym- 
bols: (100) or (010) direction). The corresponding volumes 
are also given; Vb is the volume of LaCoC>3 at room tem- 
perature and ambient pressure. Multiple points for the same 
compound are results from different structural data sets. 



HS and a LS-LS pair is 



5E Lli = (AH S - HS 



A LS-LS \ 
— ^SE J 



For C/ avg = 5.5 eV, we find 5E LR < (energy gain if 
a HS-HS is formed) for systems with larger ionic radius 
and SE 1 ^ > (energy loss) for systems with smaller ionic 
radius. The energy cost of a HS-HS pair with respect to 
two isolated HS ions is 



SE 



A£ S H E S - HS 



A£ S L 1" LS - 2A£ S H E S - LS . 



We find that <5Ase > for all systems and all parameter 
ranges. For £/ avg = 5.5 eV, it decreases from 70 to 30 meV 
along the series; for fixed Ln, it increases with decreasing 



Based on these results, we can build a model which 



describes the system. We introduce pseudospin operators 
<7 Z , and identify LS sites with spin down and HS sites 
with spin up^ The interactions between spins yield a 
Ising-like model 



H ee = J2HT)ai + lj2 J(T 

i 



% ~3 
(7 . 



— nf, and (ji) are near neighboring lattice sites. The 



Here ni — nf 1 + = 1 is the occupation per site, a\ 



parameter 



h(T) = \e^{T) + \q5E™, 



where q is the coordination number, plays the role of an 
external field. The temperature enters explicitly only in 



r h (T) = -3J aV g + 30AJ avg + 2A ljM (T) 

The coupling is J = |<5Ase > (antiferro). In static 
mean-field a LS homogeneous state is given by the solu- 
tion of the self-consistent equation 

(4) = tanh(-MT) - qJ{a\))p. 



At low temperature h(T) > is large and the system is 
in a fully polarized ferro LS state. Increasing the tem- 
perature, h(T) decreases linearly (Fig. O allowing the 
formation of some HS sites. We find that <5i? LH increases 
with decreasing ionic radius; the crystal-field is instead 
maximum around Ln=Eu (Fig. [3]). An increase of ac- 
tivation energy from 1200 K in LaCoOs to 3200 K in 
EUC0O3 has been indeed reported^ 

When HS states form, covalency-driven lattice relax- 
ation effects^ change the environment. In our model, 
if we take into account only two-sites interactions, this 
effect enhances J — » J' = J + J IC \, thus favoring a AF 
HS/LS phase, described in mean field by 

(4)=tanh(-/ l (T) + gJ'(^))/3. 

An overestimate of J re i can be obtained by using the 
difference between HS and LS in empirical ionic radii 
(friR ~ 0.06A) and the fact that A varies linearly when 
changing the Co-0 distance by a small amount; we can 
estimate the slope from our results (Fig. [3]). We find 
J rc i ~ j2SA, with 26A ~ ±2A'5r m ~ 90 meV, decreas- 
ing to 80 meV at high temperatures, i.e. comparable with 
the super-exchange term 5Ase- 

If h(T) = 0, eventually an antiferro ordered HS+LS 
lattice forms; the critical temperature for such a state can 
be estimated using mean-field theory, UbThs+ls = qJ '■ 
If J re] = 0, for (7 avg = 5.5 eV T HS+LS ~ 1200 K for 
Ln=La, decreasing to 500 K for Ln=Lu; J lc \ enhances 
Ths+ls of about 1300 K. 

The critical temperatures we obtain are very large. 
However, although h(T) — > with increasing tempera- 
ture, it likely remains comparable to J' till very high 
temperatures; thus the lattice stays disordered, and per- 
haps even close to a LS ferro solution in a large temper- 
ature range, with (a\) (i.e., the occupation of LS states) 
decreasing linearly with increasing temperature. Indeed, 
in LaCoC>3 it has been reported that the fraction of HS 
sites increasing slowly from 0.1 to 0.4 increasing the tem- 
perature from 100 to 700 K^. Hysteresis effects could arise 
because of the lattice relaxation. 

Up to here we have assumed that in HS-HS pairs Co 
HS ions stay paramagnetic and paraorbital. However, 
a specific molecular state could form as well as a spin 
singlet or triplet. We find that orbitals in the degen- 
erate ti g states play a minor role, because the energy 
differences of different orbital states with respect to the 
orbital average is small. Instead, the formation of a mag- 
netic singlet could play a role. In order to estimate the 
associated energy gain, we calculated the coupling con- 
stant Jse of the Heisenberg interaction Jse S, • Sj for 
a HS pair. We find that Jse is AFM for all systems 
(see Fig. [5]). In LaCoC-3, for L7 avg = 5.5 eV we find 
Jse ~ 7 meV at room temperature and ambient pressure; 
it increases to 10 meV if the pressure rises to 8 GPa, and 
slightly decreases with increasing temperature. Reduc- 
ing the ionic radius the anisotropy of Jse increases, but 
the average magnetic coupling decreases down to 8 meV 
for Ln=Lu. Next we assume that in all system, due to 
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FIG. 7: (Color on-line) Orbital-resolved LDA+DMFT spec- 
tral functions (ti g and e g ) of LaCoC>3 for decreasing tempera- 
ture, approaching the insulator to metal transition from above 
(T/Tim ~ 4 -> 2). Upper panels: 17 = t7 av g = 3 eV, P=0 GPa 
and 8 GPa. Lower panels: U = (7 a v g = 4.5 eV, P=0 GPa and 
8 GPa. 



the spin-orbit coupling, the effective magnetic moment 
p = \J J(J + 1) ~ \f2\ then, the magnetic energy gain as- 
sociated with the formation of a singlet is — | Jse, which 
in LaCoOa is ~ — 18 mcV. This energy gain partially re- 
duces <5Ase an d therefore J, which would change sign 
around the end of the series if U avg ~ 5.5 cV; <5Ase 
increases however rapidly with decreasing U avg . The 
smaller Ths+ls the more favorable is a disordered phase, 
reinforcing the conclusions of the paragraph above. 

However, the magnetic coupling JgE is sizable, and 
magnetic interactions would be crucial in a homogeneous 
HS-HS state. In mean-field theory, the critical temper- 
ature T N ~ | J SE J(J + 1) should then be ~ 500 K for 
LaCoC>3, decreasing in the series to 400 K and increasing 
under pressure to about 700 K. The mean-field values 
are overestimates; quantum-fluctuation effects, further 
crystal-field splittings reducing the effective magnetic 
moment, and the zero-field splitting, reduce Tn. Still, 
even taking reduction factors into account, many cobal- 
tates, if in a homogeneous HS state, should be magnetic 
because the spin-state transition occurs at Tgg < TV; in 
the presence of lonely HS-HS pairs, anti-ferromagnetic 
short range correlations could be detected. On the other 
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FIG. 8: Total LDA+DMFT spectral function (T/T IM ~ 3) 
for LaCoOs; t7 a vg = 3 eV (light) and U avg = 4.5 eV (dark). 

hand, if the number of ions thermally excited to a HS 
state is small till high temperature and the phase is disor- 
dered, weak FM short-range correlations, as reported in 
dynamic neutron-scattering experiments, could perhaps 
be triggered by a double- exchange like mechanism,— or 
even alone by the small but finite zero-field splitting 



C. High-temperature phase 

The transition observed at Tim is usually ascribed to a 
scmiconductor-to metal transition, the nature of which is 
hotly debated. There are indications that strong correla- 
tions play a crucial role, and that the spin-state transition 
is also a crossover from a spin-gapped insulator to a Mott 
insulator.— ^ In the absence of firm evidence of the forma- 
tion of super lattices, we analyze to what extent the insu- 
lator to metal transition can be understood in terms of a 
Mott transition within the homogeneous phase. Several 
works indeed suggest that the high-temperature phase 
could be a metallic homogeneous IS state.i 7 -^ To dar _ 
ify if this can be the case, we solve the Hamiltonian 
(III) by means of the LDA+DMFT approach^ we per- 
form calculations for LaCoC>3 (Figs. [3 andHJ), for which 
spectral properties are best known. The experimental 
gap is ~ 0.2 eV at low temperatures and slightly larger 
(~ 0.5 eV) above 300 K; to be consistent with experi- 
mental gaps, we vary /7 avg between 3 eV and 5 eV. 

Fig. [7] shows the LDA+DMFT results for (7 avg = 3 eV 
and the ambient pressure structure for temperatures ap- 
proaching Tim (panels (a) and (b)). We find a ho- 
mogeneous high-spin state; thus we do not find sup- 
port for the propose d 17 ' 22 homogeneous IS state at high- 
temperature. Remarkably, we find metallic t2 g and insu- 
lating e g spectral- functions, i.e., an orbital-selective Mott 
state. This can be understood from the fact the t^ g 
electrons have a larger orbital degeneracy^ although a 
smaller band- width (FigJS]), while the e g electrons are 
half-filled in the HS state, and therefore the exchange 
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FIG. 9: (Color on-line) LnCoOa series: Evolution of the e g 
(squares) and ti g (circles) bandwidth with pressure (Ln=La, 
R3c structure), temperature (filled symbols: Ln=La (R3c 
structure) and empty symbols: Ln=Y), and ionic radius Ln 3+ 
(right panel). 

coupling Jav g effectively enhances the Coulomb repul- 
sion. Decreasing the temperature towards Tim a pseudo 
gap opens in the e g spectral-function, and eventually an 
insulating regime with a small gap and a corresponding 
bad metal behavior appears. We find that the multiplct 
positions in the spectral function are close to the atomic 
limit high-spin curve in Fig. O We repeat the calcula- 
tion for the 8 GPa structure (Fig. [71 panels (c) and (d)); 
this system has a larger crystal-field A avg , very close, for 
the chosen exchange parameters, to the low spin to high- 
spin ground-state crossover in the atomic limit. Com- 
pared to the ambient pressure case, the e g gap is smaller, 
while, at fixed temperature, the ti g spectral function re- 
mains metallic with higher low-energy density of states. 
These results indicate that, as far as the state remains 
HS, at 8 GPa the insulator to metal transition should oc- 
curr at lower temperature than at ambient pressure, i.e. 
the system should stay metallic. Experiments show how- 
ever that LaCoOs returns to the insulating state under 
pressurej 24 ' 26 Our results exclude that this can happen 
in a HS scenario, and are instead compatible with the 
suggestion^ that the metal-insulator transition observed 
under pressure is driven by a spin-state transition. 

Increasing U avg to ~ 4.5 eV we find (Fig. [7J panels 
(e), (f), (g), (h)) that a real gap opens in the t2 g spec- 
tral function, and both e g and t2 g spectral functions are 
insulating even at high temperatures. The <2g gap is 
small, still compatible with the bad metal behavior ob- 
served for T > Tim- Furthermore, calculations for the 
high-pressure structure show that the t2 g spectral func- 
tion does not change much while the e g gap is reduced. 
This can be understood observing that the e g band-width 
increases substantially more that the ti g band-width in- 
creasing pressure (Fig. [5]). 

The total spectral function at ambient pressure is 
shown in Fig. [SJ We find three negative energy peaks 
(A,B,C), and two positive energy features (D,E). For 
£4v g ~ 4.5 cV peak A is at —7.5 cV, peak B at —4.7 cV, 
~ 2 eV above, and peak C at ~ —1.9 eV, features all 
observed in XPS, UPS, or photoemission dataJ£i2L£2i2i 



The spectra for U avg ~3 eV and f/ avg =4.5 eV differ in 
particular in the photoemission part; reducing J7 avg to 
3 eV, we find that the spectrum is moved almost rigidly 
toward the right by about 2 eV; the first and second fea- 
tures move to -5.5 eV and -3.5 eV respectively, while 
the lowest energy peak moves to —0.5 cV, very close 
to the Fermi edge, to partially merge with the 0.8 eV 
peak; finally, spectral weight moves from B to C and 
A. While the exact positions of peaks A, B, and C shift 
with J7 avg , the overall shape of the spectral function ap- 
pears in line with XPS, BIS, XAS and PES data at room 
temperaturej 63 i 70 ' 71 The positive energy features D and 
E at ~ 1 eV and ~ 3 eV are reflected in the form of the 
XAS and BIS spectral 

Thus the insulator to metal transition, as described 
in a homogeneous scenario, and the spin-state crossover 
exhibit different trends; the spin-state crossover is con- 
trolled by small changes in A, exchange anisotropy and 
super-exchange, parameters which change sizably with 
decreasing ionic radius. If the homogeneous HS is pop- 
ulated, a Mott insulator to bad metal transition can oc- 
cur, at a temperature which depends strongly on the ti g 
band-width and crystal-fields. 



IV. CONCLUSIONS 

We have studied the nature of the spin-state and metal- 
insulator transition in LnCo03 cobaltates. We show 
that a low-temperature intermediate spin-state scenario 
is unlikely in the full LnCoOs series. We show that the 
spin-state transition is controlled not only by the cu- 
bic crystal-field A and the average Coulomb exchange 
Javg but, surprisingly, also by the Coulomb exchange 
anisotropy AjJ avg and by super-exchange. We propose a 
simple Ising-like model to describe the spin-state transi- 
tion. We find that lattice relaxation and super-exchange 
yield anti-ferro coupling of the same order, which com- 
pete with the crystal field. Our model qualitatively ex- 
plains the trends observed in spin-state transitions in the 
LnCoOa scries. By using the LDA+DMFT approach, 
we show that in LaCoOs a high-temperature homoge- 
neous intermediate spin state is unlikely. Within a ho- 
mogeneous HS state, we find that the HS metal-insulator 
transition has a different nature than the spin-state tran- 
sition, as it mostly depends on the t2 g states and it could 
be orbital-selective. 
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